Signal Amplification in NbN Superconducting Resonators via Stochastic Resonance 
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We exploit nonlinearity in NbN superconducting stripline resonators, which is originated by local 
thermal instability, for studying stochastic resonance. As the resonators are driven into instability, 
small amplitude modulated (AM) microwave signals are amplified with the aid of injected white 
noise. The dependence of the signal amplification on the modulation amplitude and the modulation 
frequency is examined and compared with theory. 

PACS numbers: 74.40.+k, 02.50.Ey, 85.25.-j 



O ■ 

o ; 

u • 

Of 

3 ; 



-a 
c 

o 
o 



> 
in 
in 
in 
\o 
o 

O 

-i— > 
c3 



i 

C 

O 

o 



X 



The notion that certain amount of white noise can ap- 
preciably amplify small periodic modulating signals act- 
ing on bistable systems, generally known as stochastic 
resonance, has been over the last two decades of a great 
interest 0, 0, 0, 0, 13 ■ It has been applied for instance 
to account for the periodicity of ice ages occurring on 
earth by Benzi and collaborators 0, • It was demon- 
strated experimentally in several electronic devices such 
as ac-driven Schmitt triggers Q . It has a growing role in 
explaining some important neurophysiological processes 
in neuronal systems 0, IToj . In addition, it was used to 
amplify small signals in various nonlinear systems, e. g. 
the intensity of one laser mode in a bistable ring laser 
the magnetic flux in a superconducting quantum 
interference device |l2l Il3| . and even more recently, a 
small periodic drive of a nanomechanical oscillator [14j . 
In this work we exploit the metastability exhibited by 
our nonlinear NbN superconducting stripline resonators 
in order to amplify a small applied sinusoidal modulat- 
ing signal using stochastic resonance. Moreover, we apply 
stochastic resonance measurements in order to examine 
the dependence of the signal amplification generated by 
the resonator system on both the modulation frequency 
and amplitude. 

The layout of the center layer of the NbN stripline res- 
onator that was used is shown at the top-right corner of 
Fig. |2 Fabrication details as well as nonlinear charac- 
terization and modeling of these resonators can be found 
in Refs |l5lll6||. 

As it was mentioned earlier, one of the necessary con- 
ditions for the demonstration of stochastic resonance 
in a nonlinear system is the existence of metastablc 
states, which the system can hop between with the aid of 
stochastic noise and a suitable modulation drive. Thus, 
in order to set a possible working point of the resonator 
at the metastable region, two preliminary hysteresis mea- 
surements were performed. In one measurement exhib- 
ited in Fig. ^ (a), forward and backward frequency 
sweeps of the reflection parameter Sn were measured for 
the fundamental mode of the resonator at /o — 2.57 GHz. 
The bidirectional frequency sweeps form two hysteresis 
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FIG. 1: (Color online). (a) Forward and backward fre- 
quency sweeps applied to the first mode of the resonator at 
~ 2.57 GHz . The sweeps exhibit hysteresis loops at both sides 
of the resonance line shape. The plots corresponding to dif- 
ferent input powers were shifted by a vertical offset for clarity, 
(b) Reflected power hysteresis measured at a constant angu- 
lar frequency of u) p = 2n ■ 2.565 GHz which resides within the 
left-side metastable region of the resonance. For both plots 
the black (dark) line represents a forward sweep whereas the 
green (light) line represents a backward sweep. 



loops at both sides of the resonance line shape at which 
the resonator becomes bistable. In another measurement 
shown in Fig. H(b) the frequency of the pump f p was set 
to 2.565 GHz positioned at the left side of the resonance, 
while the input power was swept in the forward and the 
backward directions. A hysteresis loop of the reflection 
parameter is apparent in this measurement as well, this 
time along the power axis. Thus, the working point was 
set to f p — 2.565 GHz, Pq = —21.5 dBm, while the ap- 
plied modulation drive was a sinusoidal AM signal with 
a modulation amplitude A mo ^ = 0.27. 

A schematic diagram of the experimental setup em- 
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FIG. 2: Schematic drawing of the experimental setup used to 
measure stochastic resonance. The microwave signal genera- 
tor and the local oscillator at frequency f p were phase locked. 
The layout of the resonator is shown at the top-right corner. 



ployed to measure stochastic resonance is depicted in Fig. 
121 A coherent signal Po cos(w p t) with angular frequency 
Lu p = 2tt ■ fp is AM modulated using a sinusoidal genera- 
tor with an angular frequency f2. The modulated signal 
is combined with white noise and injected into the res- 
onator. The white noise which is generated using a noise 
source is amplified using an amplifying stage and tuned 
via an adjustable attenuator. The reflected power off 
the resonator is mixed with a local oscillator with fre- 
quency fp and measured simultaneously in the time and 
frequency domains via an oscilloscope and a spectrum 
analyzer respectively. Thus, the input signal power fed 
to the test port of the resonator (after calibrating the 
path losses) reads 



Pin(t) - Po [1 + A mod sin (Qt + 99)] cos(^pt) + £(t), (1) 

where £(t) denotes a zero- mean Gaussian white noise 
(£(*)) = 0, with autocorrelation function (£(t)£(t')) — 
2DS(t — t'), where D is the noise intensity. For conve- 
nience, we choose the phase of the periodic drive to be 
zero tp — 0. 

For small amplitudes of the modulation signal A mo< i <C 
1 and in steady state conditions, the reflected power off 
the resonator P r (t) measured via a homodyne detection 
setup at the output, can in general, be written in the 
form HQ 



Pr(t) = Po + E A n( D ) C0S + MD)) + C(t), (2) 

n=l 

where P r {t) is a superposition of a dc component Pq, 
periodic functions of time with integer multiples of the 
fundamental angular frequency fi, and £ r (i) the noise at 



the output. The coefficients A r n (D) and <j> n (D) represent 
the reflected power amplitude and the phase lag of the 
n-th harmonic as a function of the noise intensity respec- 
tively. Calculating the phase-averaged spectral density of 
P r (t) yields 



S(w) = 27rP r f5(w) (3) 

OO 

+ 7T A n ( D ) I 6 ( w -nQ) + S(cj + nfl)] + S N (w), 

n=l 

which is composed of a delta spike at dc (uj = 0), delta 
spikes with amplitudes A r n (D) centered at uj — ±nf2, 
n = 1,2,3..., and a background spectral density of the 
noise denoted by Sn(u>). 

As it is known, one of the distinguished fingerprints 
of stochastic resonance phenomenon is a peak observed 
in the signal to noise ratio (SNR) curve as a function 
of the input noise intensity D, corresponding to some 
nonzero intensity -Dsr- This counterintuitive amplifica- 
tion in SNR curve is generally explained in terms of co- 
herent interaction between the modulating signal and the 
stochastic noise entering the system. 

In this framework a SNR for the n-th harmonic can be 
denned as 



7lQ + AiU 



SNR„ = 2 



lim 



S (uj) did /Sn (nfl) 



ifi— Aa 



^Al(D) 
Sn (nfi) 

where SNR of the fundamental harmonic corresponds to 
n = 1. 

In Fig. |3 three SNR curves as a function of noise 
intensity are shown, corresponding to the fundamental 
n = 1 and the odd harmonics n = 3 and n = 5. All three 
curves display a synchronized peak in the SNR around 
P>sr of about 2 fW/Hz. 

Typical results of stochastic resonance measured in the 
time and the frequency domains are shown at the left and 
right sides of Fig. ^respectively. Panels (a) and (b) cor- 
respond to low noise levels below the critical value. Panel 
(a) shows the reflected sinusoidal at /n = 1 kHz without 
jumps. Panel (b) shows the reflected sinusoidal contain- 
ing a few arbitrary jumps. Whereas, panel (c) which 
corresponds to the critical stochastic noise Dsr exhibits 
one jump in the reflected signal at every half cycle. That 
is, the time-scale matching condition for stochastic reso- 
nance given by t(Dsr) = Th/2, is fulfilled for this noise 
intensity, where Tq = 27r/f2, and t(D) is the metastable 
state lifetime corresponding to the noise intensity D. A 
possible theoretical linkage between the transition rate 
1 /t exhibited by these resonators and their physical non- 
linear mechanism will be established in a future publica- 
tion [l^. In panels (d) and (e) on the other hand, the 
case of noise levels higher than Dsn are shown. In panel 
(d), the coherence between the modulating drive and the 
noise is lost, and multiple jumps are observed during each 
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cycle. Finally, in panel (e) , the high noise almost screens 
the signal and the induced jumps are totally uncorre- 
cted. It is also apparent that the reflected modulated 
signals shown in panels (c) and (d), in addition to be- 
ing amplified, exhibit a rather rectangular shape. Such 
distortion to the sinusoidal shape at the output may orig- 
inate from a dispersive character of the system response 
as was suggested in Refs. [a, Ell- 
in the frequency domain displayed at the right side 
of panels (a)-(e), the fundamental and the harmonics up 
to n — 5 are shown, as the noise intensity is increased. 
The appearance of even harmonics in the frequency spec- 
trum indicates that the nonlinear system is asymmetric, 
in contrast to symmetric potential systems where only 
the odd harmonics emerge 0, ^| . 

Furthermore, if one defines a spectral amplification of 
the ra-th harmonic spectral density relative to that of the 
fundamental measured for the noiseless case, one gets 



FIG. 3: SNR curves measured for the fundamental n — 1 and 
the odd harmonics n — 3, n — 5, as a function of the input 
noise intensity. A peak in the SNR is detected at around 
Dsn of 2 fW / Hz. 
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FIG. 4: (Color online). Panels (a)-(e) exhibit typical snap- 
shots of the reflected modulated signal in the time domain 
(left side) and in the frequency domain (right side) as the 
input noise intensity D is increased. Panels (a)-(b) corre- 
spond to noise intensities below Dsr- Panel (c) corresponds 
to a noise intensity of Dsr. = 0.94 fW/Hz. Panels (d)-(e) 
correspond to noise intensities higher than Dsr- The dotted 
sinusoidal line drawn in the time domain represents the mod- 
ulation signal. 



Vn (D)=A r n (D)/AI(D = 0). 



(5) 



where the amplitude A\ (D — 0) is assumed to be propor- 
tional to the input power of the modulation signal up to 
some constant reflection factor. Since the different har- 
monics of the signal are strongly correlated, we choose in 
the following figures to exhibit the amplification of the 
third harmonic instead of the fundamental one n = 1. 
This is mainly because the first harmonic resides at the 
frequency of the modulation source signal and part of it 
bypasses the resonator. 

In Fig. the amplification parameter 773 of the third 
harmonic as a function of the input noise intensity was 
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FIG. 5: Third harmonic amplification 773 as defined in Eq. 
plotted as a function of input noise intensity. The differ- 
ent plots correspond to fn = 5, 6, 7, 8, 9 kHz modulation 
frequencies. The inset shows 773 gain as a function of fn, cor- 
responding to a constant input noise intensity of 1 fW / Hz . 



4 



0.025 




2 4 6 8 10 12 14 16 

Noise intensity [fW/Hz] 



FIG. 6: Third harmonic amplification 7/3 as defined in Eq. 
plotted as a function of input noise intensity at fn — 1 kHz. 
The different plots correspond to an increasing modulation 
amplitude A mod = 0.09, 0.18, 0.27, 0.35. 



measured at a fixed modulation amplitude but for dif- 
ferent AM drive frequencies /n = 5, 6, 7, 8, 9 kHz. As 
it is apparent from the figure 773, amplification is non- 
monotonic upon increasing the forcing frequency /o . The 
amplification graphs rise up with frequency up to 7 kHz 
and afterwards decay down. A cross section of 773 param- 
eter at a constant noise intensity of 1 fW / Hz is shown 
in the inset as a function of fa varying in the range 
0.5 — 10 kHz . A peak in 773 is observed at fn = 7.5 kHz . 
Similar amplification behavior was observed for the fun- 
damental and for the fifth harmonic as well. However, 
this non-monotonic behavior of 7/3 as a function of the 



forcing frequency somehow departs from the monotonic 
decrease presented for example in Refs. 0,0. This may 
be due to the nature of the nonlinearity governing our 
system |l6[ which can not be explained in terms of Duff- 
ing oscillator nonlinearity, that was assumed in the cited 
papers. A theoretical prediction of a non-monotonic SNR 
versus fn was introduced though in Ref . 20] , where the 
case of two square potential wells separated by a square 
barrier was considered. 

In contrast, when the modulation amplitude is in- 
creased at a fixed forcing frequency, the resultant ampli- 
fication is expected to decrease, mainly due to increased 
nonlinear effects at higher modulation amplitudes, which 
tend to suppress stochastic resonance phenomenon Q. 
A good confirmation to this hypothesis is presented in 
Fig. [f)| which displays decaying amplification curves, as a 
function of the noise intensity, upon increasing the modu- 
lation amplitudes A mo a = 0.09, 0.18, 0.27, 0.35 at a fixed 
forcing frequency of /q = 1 kHz . 

In conclusion, nonlinear NbN superconducting res- 
onators have been shown to exhibit stochastic resonance 
when driven into the bistable region. Amplification of a 
slowly varying AM signal carried by a microwave pump 
can be achieved by establishing a resonant cooperation 
between the modulating signal and the injected stochas- 
tic noise. Such amplification scheme may be useful in the 
detection and amplification of small AM signals in com- 
munication area. Moreover, stochastic resonance phe- 
nomenon can be employed to some extent, to "probe" 
some important characteristics of the nonlinearity mech- 
anism of the resonators 17]. 
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